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Abstract SARAS is a correlation spectrometer purpose designed for preci- 
sion measurements of the cosmic radio background and faint features in the 
sky spectrum at long wavelengths that arise from rcdshifted 21-cm from gas in 
the reionization epoch. SARAS operates in the octave band 87.5-175 MHz. We 
present herein the system design arguing for a complex correlation spectrom- 
eter concept. The SARAS design concept provides a differential measurement 
between the antenna temperature and that of an internal reference termina- 
tion, with measurements in switched system states allowing for cancellation 
of additive contaminants from a large part of the signal flow path including 
the digital spectrometer. A switched noise injection scheme provides absolute 
spectral calibration. Additionally, we argue for an electrically small frequency- 
independent antenna over an absorber ground. Various critical design features 
that aid in avoidance of systematics and in providing calibration products for 
the parametrization of other unavoidable systematics are described and the ra- 
tionale discussed. The signal flow and processing is analyzed and the response 
to noise temperatures of the antenna, reference termination and amplifiers is 
computed. Multi-path propagation arising from internal reflections are consid- 
ered in the analysis, which includes a harmonic series of internal reflections. 
We opine that the SARAS design concept is advantageous for precision mea- 
surement of the absolute cosmic radio background spectrum; therefore, the 
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design features and analysis methods presented here are expected to serve as 
a basis for implementations tailored to measurements of a multiplicity of fea- 
tures in the background sky at long wavelengths, which may arise from events 
in the dark ages and subsequent reionization era. 

Keywords Astronomical instrumentation, methods and techniques • 
Methods: observational • Cosmic background radiation • Cosmology: 
observations • Dark ages, reionization, first stars • Radio continuum: ISM 

1 Introduction 

Following recombination at redshift z ~ 1100, the propagation of the cosmic 
microwave radiation (CMB) through the 'dark ages' and 'cosmic dawn' results 
in deviations of its spectral form from a Planck blackbody spectrum. The early 
signatures arise owing to radiative transfer interactions of the CMB with the 
spin-flip 21-cm transition in the neutral hydrogen content of the universe, 
which depends on the spin temperature evolution of the gas. Subsequently, 
in the post-reionization era, free-free emission from the ionized intergalactic 
medium as well as Compton scattering of CMB photons by the hot electrons 
in this gas imprint further deviations from the Planck form. All of the above 
traces, however, constitute a small fraction of the total cosmic radio back- 
ground at long wavelengths: the relatively more recent history of star forma- 
tion and AGN activity in galaxies cumulatively manifest as the radio source 
counts and constitute the bulk of the extragalactic radio background. 

Cosmological expansion and the consequent redshift results in a movement 
of the cosmological radiation signatures to longer wavelengths; therefore, pre- 
cision measurements of the spectrum of the cosmic radio background (CRB) 
at long wavelengths is an important probe in modern cosmology. 

Measurements to date suggest that the extragalactic radio background is 
the CMB plus a power-law spectrum that has a brightness temperature of 
(1.2 ± 0.09) x {v/l GHz)- 2 - 60±0 04 K (Fixsen et al. 2009). Surprisingly, the 
sky brightness corresponding to discrete radio sources detected in the deepest 
surveys to date account for only a fraction of the extragalactic radio back- 
ground, even after excluding the CMB (Seiffert et al. 2011; Singal et al. 2010). 
Improved measurements of the radio background and, in particular, the spec- 
trum at long wavelengths where errors are relatively larger, are important in 
estimating the spectrum of the unexplained part and thereby constraining the 
sources of this cosmic radiation. 

Free-free emission from ionized gas in the post-reionization era adds to 
the radiation background, and is usually characterized by the optical depth 
to free-free emission Yg (Bartlett and Stebbins 1991). Energy release at early 
times — for example, as a consequence of radiative decay or annihilation of 
relic particles — would be expected to manifest as distortions to a Planckian 
spectrum of the /x form (Sunyaev and Zeldovich 1970). Constraints to Yg 
are derived by fits of models for the evolution in the gas to measurements 
of the radio background; constraints on [i are derived from fits of spectra 
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with a Bose-Einstein photon occupation number where the deviation from 
Planckian form is parametrized by /j,. Deviations in the background brightness 
temperature at frequencies below about 1 GHz best constrain these terms 
and, therefore, provide a motivation for improved measurements of the long 
wavelength spectra of the extragalactic radio background. 

Perhaps the key scientific motivation today for any measurement of the 
spectrum of the radio background at long wavelengths is the cosmological 
epoch of reionization. The evolution of the spin temperature of the gas, as the 
coupling between hyperfine level populations with gas kinetic temperature and 
radiation background changes with cosmic epoch, is expected to manifest as 
a distinctive signature in the radiation spectrum (Pritchard and Loeb 2008). 
Significant deviations in the spin temperature of the gas from the tempera- 
ture of the ambient radiation background commence early, at redshift beyond 
100, as the spin temperature is driven below the radiation temperature by 
the gas cooled by cosmological expansion. However, this coupling effectively 
disappears by about redshift 30. Spectral features associated with the forma- 
tion of the first condensed objects are expected to appear at redshift about 
20, when the coupling of spin to the gas kinetic temperature is renewed by 
the Wouthuysen-Field coupling (Wouthuysen 1952; Field 1958) induced by the 
Lyman-a from the first light. Subsequent heating of the gas and reionization 
result in a transformation of the gas visibility from absorption to emission, 
and a diminishing of the gas emissivity in redshifted 21-cm as the neutral 
fraction vanishes; this final step in the reionization epoch is expected to occur 
progressively between redshifts of about 13 and 6 (Zahn et al. 2012). 

Reionization is expected to be effectively complete by redshift of 6, cor- 
responding to an observing frequency of about 200 MHz, and hence all these 
features associated with the cosmological evolution of the gas is at lower fre- 
quencies. However, the magnitude of the signal corresponding to the epoch 
of reionization is not expected to exceed 0.1 K and is at least three orders of 
magnitude smaller than the foreground. The astrophysical evolution including 
the formation of the first compact objects, stars, galaxies, and the ionizing 
radiation from these sources are poorly constrained: therefore, precision mea- 
surements of deviations in the background spectrum from Planckian form are 
means to improved modeling of this era. 

We have evolved a design for a correlation spectrometer for the measure- 
ment of the absolute spectrum of the cosmic radio background in the frequency 
range 87.5-175 MHz. The system has been given the acronym SARAS, for 
Shaped Antenna measurement of the background RAdio Spectrum. SARAS 
deployed in the Gauribidanur Observatory about 80 km north of Bangalore 
in India is shown in Fig. [T] In this paper we present the design of SARAS, 
with a focus on features — design strategies and calibration methods — that are 
novel and relevant for wideband measurements of background spectra. The 
work presented herein represents progress towards building a system capable 
of a useful measurement of spectral signatures of the EoR in the cosmic radio 
background. 
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Fig. 1 The SARAS Fat-dipole correlation spectrometer 

We discuss the SARAS system configuration in Section 2 including details 
of the antenna as well as the analog and digital receivers. Because the work 
presented herein is a pioneering effort in that it presents a complex correlation 
approach to solving for total power spectra, we have chosen to present the 
performance in stages beginning with a simplified treatment and then going 
on to develop the understanding considering non-ideal behavior of increas- 
ing complexity. In Section 3 the measurement equations, signal processing 
and calibration strategies are first described assuming an ideal system perfor- 
mance. Modified measurement equations taking into consideration first-order 
non- idealities are discussed in Section 4.1, with a more complete treatment 
considering higher order corrections presented subsequently in Section 4.2. 
The approach to modeling of SARAS data is discussed in Section 5. 



2 SARAS system configuration 

We present here a system level design for a correlation spectrometer based 
receiver for the measurement of the absolute spectrum of the sky brightness. 
The signal from the antenna is split early on, amplified by identical receiver 
chains, and processed by a digital spectrometer. The two signals are sampled, 
the sampled data are discrete Fourier transformed, and data in corresponding 
frequency channels multiplied and the complex correlations integrated to form 
cross spectra over the observing band. 

The primary motivation for adopting a correlation receiver configuration 
is that it makes possible switching to cancel spurious spectral features arising 
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from additive correlator-based errors. Digital receivers, which include digital 
signal processing devices that perform complex algorithms involve signals that 
are clocked over a range of rates and generate a broad spectrum of unwanted 
emissions that propagate across transmission paths via radiative and conduc- 
tive coupling. The high-speed clocking of digital signals also generates noise in 
system electronics ground. Common mode coupling of these signals into the 
analog signal paths leading to the samplers results in spurious correlations, 
which may have frequency structure, and need to be canceled. The system de- 
scribed here switches the antenna signal so that the sky signal appears in the 
complex spectrum with a sign that flips on alternate switch positions, where 
as the spurious additive correlations would be expected to remain unchanged, 
making it possible to subtract every set of spectra acquired in adjacent time 
intervals and in alternate switch positions and thereby cancel the spurious 
correlations. 

In a correlation receiver the antenna signal may be split prior to encoun- 
tering any amplifier and, therefore, the system response does not contain the 
receiver noise temperature. This source of noise in the signal path only appears 
in the measurement response of a correlation receiver because of unwanted 
leakage (non ideal performance) in the power splitter. However, in the SARAS 
design, the switching that cancels additive correlator-based errors also cancels 
such leakage in the signal path. 

When calibrated for the complex bandpass, the sky power arriving at the 
antenna would appear wholly in the real component of the complex spectra 
measured by the correlation spectrometer. However, reflections internal to the 
receiver chain cause signals to arrive at the samplers via multi-path propaga- 
tion. As discussed below, this is a route whereby receiver noise does appear 
in the measurement response. The multi-path propagation and asymmetries 
in the system result in the measurement manifesting complex correlations in 
which the real and imaginary components share many descriptive parameters. 
This provides an opportunity for using the imaginary components of the mea- 
sured spectra, in which no sky signals are present, to estimate parameter val- 
ues that characterize the multi-path related spectral features that contaminate 
the real spectra. These parameters, which are derived from the measurements 
themselves, may then be used in a joint modeling of (a) the components that 
constitute the system temperature, and (b) the propagation paths internal to 
the receiver, and thereby reveal the characteristics of the sky power arriving 
at the antenna. This is a secondary advantage of adopting a correlation re- 
ceiver configuration that includes a switch, and provides additional methods 
for modeling the sky power and distinguishing features owing to the cosmic 
radio background from multiplicative and additive features associated with 
internal reflections within the signal path. 

The configuration of the SARAS spectrometer is shown in Fig. [2] In the 
sections below we describe each of the components in this signal path and 
discuss the rational in the design. 
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Fig. 2 The configuration of the SARAS spectrometer 



2.1 The Antenna: a frequency independent fat-dipole 

The antenna samples the sky brightness and couples the sky power to the 
signal path. The sampling is defined by the response to the distribution of sky 
brightness and by the frequency dependence of this response. The coupling 
efficiency of an antenna is related to the ohmic loss and the impedance match 
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that the antenna offers to the electromagnetic wave arriving in free space. The 
sky brightness has structure on a range of scales and frequency dependence 
in the antenna pattern will result in a spectral response that varies with fre- 
quency, which will depend on the frequency dependence of the antenna and 
the variation in sky brightness in the direction of the antenna beam. For this 
reason a frequency independent antenna is desirable. 

An antenna that has physical size large compared to the wavelength of 
operation is also undesirable since any departures from frequency independent 
performance owing to manufacturing tolerances will result in spurious higher 
order components in the spectrum (owing to brightness structure on the sky) 
that may confuse cosmological signatures. Antenna structures that are electri- 
cally long — like a frequency independent spiral antenna, for example — would 
have return loss and, therefore, a transfer function that would have high order 
structure. High order frequency structure in the antenna transfer function will 
lead to spurious confusing spectral structure if there are errors in the calibra- 
tion of the antenna response; therefore, it is desirable that the antenna return 
loss is smooth so that the system is tolerant to antenna calibration errors. 
The above considerations require that the antenna be frequency independent, 
of small overall electrical size, and that the conductive parts be electrically 
small. For the same reasons, it is desirable that any associated balun also be 
electrically small. 

We use a frequency independent antenna that is electrically small: the an- 
tenna adopted for the SARAS spectral radiometer is a variant of the profiled 
fat-dipole (Raghunathan et al. 2012) designed to operate in the 87.5-175 MHz 
band. The dipole elements have sine profiles and the antenna has been mea- 
sured to have a beam pattern that deviates by less than 3% from a cosine- 
squared radiation pattern, which is expected for a short dipole, over the octave 
bandwidth. The arms are of square cross section and constructed using alu- 
minum sheets. In place of the coaxial-choke balun adopted by Raghunathan 
ct al. we have deployed a balun consisting of a transformer and a 180° hybrid, 
which are placed in between the arms and this module serves to transform the 
balanced out-of-phase signals from the two arms to an unbalanced signal on a 
coaxial cable. This variation is motivated by the necessity to maintain short 
transmission lines between the antenna and receiver. 

The antenna has a styrofoam base for structural stability and is placed on 
ground covered by Panashield SFA-type flat absorber tiles. The height of the 
center of the dipole is 0.81 m from the tiles, which are of Nickel-Zinc Ferrite 
composition and have a reflectivity lower than —30 dB over the band. The 
absorber ground affects the antenna temperature in two ways: (a) there is 
multi-path propagation of sky radiation to the antenna and from every sky 
direction the direct ray adds with a ray that arrives at the antenna after reflec- 
tion from the absorber ground, the vector sum results in either constructive 
or destructive interference depending on the relative path delay and phase 
change on reflection; (b) the absorber ground emits with a brightness that 
depends on the ambient temperature and the reflectivity Both these are in- 
deed frequency dependent because the reflectivity is frequency dependent. Our 
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analysis of the consequent limitations to performance shows that the emission 
from the ground is an additive contribution to the antenna temperature that is 
expected to have a broad variation over the octave band with magnitude about 
150 mK at the ends and significantly smaller at the band center, where the re- 
flectivity is a minimum. The multi-path propagation is of greater consequence 
and results in at most 1% distortion to the sky spectrum. 

The small dependence of the antenna pattern on frequency over the oc- 
tave operating band, the electrically short distance between the antenna and 
ground plane, along with the small reflectivity of the absorber tiles covering the 
ground, are together primarily a frequency dependent multiplicative gain that 
is subsumed into the calibration of the antenna gain: calibration is discussed 
below in Section 5. 



2.2 Analog receiver configuration 

The analog signal processing is in part implemented in a set of modules in 
a receiver unit at the antenna base and the remainder is in a base station 
about 100 m away where the digital receiver is also located. This separation is 
necessary to prevent any self-generated RFI from the digital system entering 
the signal path via the antenna. 

The separation of the antenna base electronics from the digital receiver and 
the analog transmission of signals across this distance could result in multi- 
path propagation of signals owing to reflections. Multi-path propagation that 
is confined to within each of the arms of the correlation receiver would effec- 
tively be a modulation of the individual band pass responses of the two arms, 
which would be expected to be calibrated as part of the complex bandpass 
calibration. However, multi-path propagation that couples signals from one 
arm to the other would result in complex correlations that would leave resid- 
ual modulations in the real and imaginary components of calibrated spectra. 
Moreover, if signals propagate along the 100-m transmission line in a reverse 
direction from the base station to the antenna and if they then couple across 
to the second arm of the correlation receiver, the spurious modulation in the 
calibrated spectra could have fine spectral structure with a period of about 
1 MHz (depending on the speed of signal transmission in the 100-m line). 
These unwanted couplings between the arms that occur at the antenna elec- 
tronics may arise because of reflections associated with impedance mismatch 
at the antenna terminals. Alternately, unwanted coupling may also arise as 
a result of radiative leakage of signal power from the transmission line that 
is picked up by the antenna. We have avoided these issues by implementing 
transmission across the 100 m distance on a pair of optical fibers, which pro- 
vides excellent reverse isolation and also eliminates the potential for radiative 
couplings between the transmission lines and to the antenna. 

As shown in Fig. [2j the signal from the antenna and a calibrated broad- 
band noise source (referred to as CAL) are fed to the direct and coupled port 
respectively of a directional coupler. The broadband noise source provides an 
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excess calibration noise when it is 'ON'. The output of the directional coupler, 
which consists of antenna noise plus calibration signal (when the CAL in ON) 
is fed to port Jl of a cross-over switch. The second input port J4 of this switch 
is fed power from a second noise source that serves as a reference (referred to as 
REF). As discussed below, the recorded spectra in the SARAS configuration 
are always difference spectra: what is measured is the difference between the 
powers from the antenna (plus CAL power when the calibration noise source 
is ON) and reference (the value of REF power depends on whether or not the 
reference noise source is on). The reference is an ambient temperature noise 
source when the second REF noise source is in the OFF state, and provides 
an excess reference power when it is in the ON state. The difference between 
spectra recorded in the two states of the reference noise source is used to de- 
rive additional constraints while modeling the systematics. The output ports 
J2 and J3 of the cross-over switch are connected respectively to the sum and 
difference ports of a 4-port power splitter, whose output ports Pi and P2 arc 
connected to independent low noise amplifiers at the head of a pair of separate 
analog receiver chains. The cross-over switch is electro- mechanical, with low 
loss and high isolation. At the antenna base, the pair of analog chains include 
a low noise amplifier and a low pass filter, which are followed by modular 
fiber-optic transmitters that include transconductance amplifiers and couple 
the electrical signal into single-mode 1550 nm optical fibers. 

The calibration noise source is a 25 dB excess noise source followed by a 
3 dB attenuator that connects the calibration noise to the coupled port of 
the directional coupler. The reference load consists of another 25 dB excess 
noise source followed immediately by a 23 dB attenuator. In both cases the 
equivalent temperatures added by the calibration noise and reference noise are 
similar to the system temperature when both noise sources are off. 

At the base station, the analog receiver includes an optical receiver followed 
by a cascade of low and high pass filters with amplifiers and attenuators in 
between. The filters limit the band between 87.5 to 175 MHz; for initial field 
testing in sites where FM transmission is a potential source of RFI, high pass 
filters limit the lower end of the band to 110 MHz. 



2.3 Digital receiver configuration 

The analog signals are sampled at 175 MHz. The sampled data streams are 
discrete Fourier transformed to decompose the signals over 1024 frequency 
channels. Complex data streams from the two receiver chains at corresponding 
frequencies are multiplied and the complex products accumulated to represent 
cross-power spectra. Autocorrelation spectra of the digital signals in the two 
receiver chains are also computed. These are all read into a laptop PC and 
stored for offline processing. The digital signal processing is performed on a 
single Virtex-5 FPGA chip. 

The digital receiver includes a GPS disciplined frequency reference that 
provides sufficient long term stability, a frequency synthesizer that generates 
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the sampling clock and a digital i/o card for generating signals that switch the 
calibration noise, reference noise and cross-over switch. All switching signals 
are sent to the antenna base via optically isolated copper transmission lines. 

The digital systems are in a shielded enclosure that is sealed from the 
environment. Cooling is achieved using heat pipes which exchange heat within 
the enclosure and to external air by forced convection using fans. The d.c. fans 
are brushless; however, in order to drive the fan rotor its position is sensed by 
a Hall-effect sensor and the stator fields are flipped by switching the current 
direction. This produces measurable broad band RFI from the digital receiver 
because a pair of fans are necessarily mounted outside the sealed enclosure: a 
noise power of order 10 K enters the antenna if such fans are located 1 m from 
the antenna without shielding. We have suppressed this contribution to be less 
than 1 mK by including RFI shielding in the air vents of the fan enclosure and 
locating the digital receiver 100-m from the antenna. 



3 Signal processing in the SARAS radiometer 

Let T a represent the antenna temperature at the mainline input port of the 
directional coupler, T ca i represents the excess calibration noise (CAL) at the 
coupled port of the directional coupler when the CAL is switched on, T re /o 
represents the reference temperature at the second port of the cross-over switch 
when the reference noise source (REF) is in 'off' state and T re t\ is the excess 
noise when REF is in 'on' state. All temperatures are referred to the mainline 
input port of the directional coupler, where the antenna would be connected. 
gsi and g^2 represent the gains of the paths from the £ port of the power 
splitter to the two samplers in the digital cross correlation spectrometer, g^i 
and gA2 are the gains of the paths from the A port to the two samplers. 
All of the temperatures and gains are functions of frequency; the <?-terms 
represent voltage gains. For clarity, we herein choose to represent the voltage 
gains between the A and £ ports and PI and P2 ports of the power splitter 
by terms that are ideally unity, and explicitly include a negative sign for the 
voltage transfer function between A and P2. These symbols, the conventions 
we adopt for denoting switch positions and states of the noise sources are 
depicted in Fig. [3j 

The measurement and calibration sequence consists of alternately switching 
the mainline output of the directional coupler and the reference noise source 
between the £ and A ports of the power splitter by means of the cross-over 
switch. At each of the positions of the cross-over switch, the calibration noise 
CAL is switched on for a time r ca ; while the reference noise source is off, 
the reference noise REF is switched on for a time r re f while the calibration 
noise is off, and for a time r src both sources of noise are off. We denote the 
'direct' switch position with the antenna connected to A port as '0', and 
the 'crossed' switch position with the antenna connected to £ port as T'. 
Complex spectra are acquired with the system cycling through a set of six 
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paths. 
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Table 1 The observing schedule with SARAS. The system cycles through the states listed 
in the rows of the table and time averaged complex spectra are recorded separately and 
sequentially corresponding to each of these states. 
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states listed sequentially in the successive rows of Table [T] r ca i , r src and r re f 
are all nominally the same time and the integration time is 0.7 s in each state. 

In the remainder of this section, we present an analysis of the signal pro- 
cessing and calibration assuming ideal impedance matching between compo- 
nents and transmission lines. In later sections we deal with the issues related 
to impedance mismatches. 

3.1 Measurement equations 

The antenna signal (or the antenna signal plus the calibration noise) and 
reference noise are connected to the sum (E) and difference (A) ports of the 
power splitter. Therefore, the cross power spectrum measured by SARAS is at 
all times the difference between the antenna temperature and reference. The 
receiver noise does not appear in the response because (a) the noise in the 
receivers in the two signal paths following the power splitter are uncorrelated 
and (b) under the assumption that the interfaces are impedance matched 
and the antenna has zero return loss the receiver noise from one arm of the 
correlation receiver does not leak into the other arm. 

The measured complex cross powers in the two switch states when both 
the reference and calibration noise sources are off may be written as: 

Pooff — —gAl9A2T a + 9El9S2T re fO + Pcor (1) 

and 

Ploff = gSigmTa — gAl9A2T r efO + Pcor- (2) 

P CO r denotes the spurious additive complex cross power spectrum that may 
be present in the measured spectrum because of undesired coupling of power 
between the transmission paths and samplers in the two arms. P cor is expected 
to be a constant and independent of the state of the cross-over switch and the 
state of the calibration and reference noise sources. The measured complex 
cross powers in the two switch states when the calibration noise alone is on 
may be written as: 



Pocal = —gAl9A2{T a + T ca i) + gsigS2T re fO + Pcor 



(3) 
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and 

Plcal = 9Sl9S2(T a + T ca i) - 9 A\9 AlTrefV + Pcor- (4) 

The measured complex cross powers when the reference noise source alone is 
on may be written as: 

Poref = -9Al9A2Ta + 9Sl9S2T r efl + Pcor (5) 

and 

Plref = 9SigS2T a - 9Al9A2T r efl + Pcor- (6) 

The subscripts '0' and '1' for the powers denote spectra for the two positions 
of the cross over switch. All of the above powers are functions of frequency. 



3.2 Cancellation of spurious additive component of the response 

In the two switch positions the antenna (and the calibration signal) and ref- 
erence are connected to the £ and A ports of the power splitter alternately. 
Therefore, in the two switch positions, the difference spectra — difference be- 
tween the signals at the £ and A ports of the power splitter — appear in the 
measurement with opposite sign. However, in both switch positions the mag- 
nitude and sign of the spurious additive contribution P cor that arises because 
of unwanted couplings of noise across the two samplers or between the pair 
of signal paths from the power splitter to the samplers remains unchanged. 
Differencing the '0' and '1' spectra corresponding to the two switch positions, 
therefore, cancels any such spurious additive contribution to the measured 
spectra. 

The difference cross power spectra are 
Poff = Pioff ~ I':, J f 

= {gSl9£2 + 9Al9A2){T a ~ T ref0 ), (7) 



Peal — Plcal Pocal 

= {9E19E2 + 9Al9A2){T a + T ca i - T ref0 ) (8) 

and 



Pref Plref P()ref 

= {9S19S2 + 9Al9A2){Ta - T re fl)- 



(9) 
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3.3 Calibration for the complex instrument bandpass 

The calibration noise source is coupled to the signal path from the antenna 
immediately following the balun. Along with the antenna signal, the calibra- 
tion noise is also connected alternately to the E and A ports of the power 
splitter, and in each switch position the complex cross-correlation spectrum is 
recorded separately with the noise source switched on and off. Neglecting the 
asymmetries in the switch and differences in the reflection and transmission 
properties of the different switch paths, the calibration noise may be used to 
calibrate the antenna temperature for the bandpass of the signal path down- 
stream of the cross-over switch. We assume here that the excess noise at the 
reference port of the cross-over switch has a flat temperature spectrum. 

Subtraction of the P f / spectrum from the P ca i spectrum gives the instru- 
ment response to the calibration noise: 

Peal - Poff = (gSl9E2 + 9Al9A2)T C al- (10) 

Since the calibration noise samples both the E and A ports of the power 
splitter, the above power spectrum that represents the excess calibration noise 
includes the mean multiplicative gains of both ports. 

The difference spectrum P ff is bandpass calibrated using {P ca i — Poff) 
and scaled by the calibration noise temperature T ca i to derive a calibrated 
spectrum in antenna temperature units. This calibration yields the desired 
antenna temperature minus the temperature of the reference port: 

T a - T ref0 = — — r x T ca i. (f 1) 

{^cal - t^off) 

It may be noted here that the gains gsi, gs2, 9ai, and g^2 cancel in this 
process. It is to make this happen that the SARAS configuration is designed 
so that the calibration noise is switched on and off in both switch states and 
the cross-over switch directs both antenna noise and calibration noise temper- 
atures via both the E and A ports of the power splitter. 

Following calibration for the complex bandpass, the real component of the 
calibrated cross spectrum represents the antenna temperature minus noise 
temperature of the reference. At frequencies where the antenna temperature 
exceeds the reference temperature the real component will be positive, where as 
over frequencies where the antenna temperature is less than the reference noise 
temperature the real component will be negative. The imaginary component 
is expected to vanish in the ideal case considered in this section. 

If the impedance match along the signal paths were excellent (ideal) then 
the reference noise is unnecessary. However, as discussed below, this second 
noise source has been included in the SARAS configuration to serve as a 
method of deriving a template for the instrument response to the reference 
load. The SARAS measurement of the sky spectrum is a differential measure- 
ment in that it is the difference between antenna temperature and that of 
the reference load; the calibration noise serves to provide a template for the 
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instrument response to the antenna signal where as the reference noise source 
serves to provide a template for the instrument response to the reference load. 
It may also be mentioned here that an additional source of significant noise 
temperature in the receiver is the noise temperature of the first low noise am- 
plifiers in the two arms of the correlation spectrometer. Although in an ideal 
system this noise would not contribute to the measured spectrum; however, 
impedance mismatches and leakage terms in the power splitter result in a re- 
sponse to receiver noise and it would be useful to have low noise amplifiers 
whose noise temperatures may be switched between nominal and high states! 
We are not aware of how such a feature may be implemented in practice. 

If we assume that the calibration noise temperature is a constant over the 
octave observing band, then T ca i is independent of frequency. The value of 
T ca i is determined by replacing the antenna with terminations at known tem- 
peratures and comparing the step in measured power with the step measured 
when the calibration noise is switched between on and off states. 



3.4 Calibration for the antenna gain 

If the antenna were ideal in that the antenna impedance matched free space 
and the ohmic losses are negligible, the antenna temperature would be the 
mean brightness temperature over the radiation pattern, where the weighting 
is the normalized radiation pattern. The frequency dependance of the gain loss 
due to impedance mismatch and ohmic loss may be calibrated by modeling 
the difference spectra when the sky brightness in the main lobe of the beam 
changes significantly, assuming that we know the difference spectrum of the 
sky brightness. We examine here considerations for the modeling of the sky 
brightness that are needed for calibration of the antenna gain. 

The brightness temperature of the sky at the center of the SARAS band is 
about 165 K in the coldest regions and the brightness rises to several thousand 
Kelvin towards the Galactic Center. The cosmic radio background towards 
most sky directions are very nearly power law spectra, with a temperature 
index that has a mean of about —2.55 (Kogut et al. 2011; Seiffert et al. 2011). 
The index varies somewhat with direction. We adopt the 408-MHz sky survey 
of Haslam et al. (1982) as representing the sky brightness distribution, and 
assume that the spectral index towards any sky direction is given by the 2- 
point spectral index computed between this image and the 150-MHz all-sky 
map synthesized by Landecker and Wielebinski (1970). 

The SARAS fat-dipole antenna has a dipole radiation pattern that aver- 
ages over a beam solid angle covering two-thirds of the sky. Adopting the above 
model for the brightness distribution over the sky and across the frequency 
band, the SARAS dipole has an antenna temperature at the band center fre- 
quency of 131.25 MHz whose variation with LST is as given in Fig. |4j The 
ground has been assumed to be a perfect absorber with brightness 300 K, 
and the variation in antenna temperature is displayed for an EW as well as 
NS orientation for the dipole. A location at +13° latitude has been assumed. 
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Fig. 4 Expected variation in antenna temperature over LST at the band center frequency of 
131.25 MHz. The red curve is for EW dipole orientation and blue curve is for NS orientation. 

When the system is used to detect spectral signatures in cosmological radio 
backgrounds such as the all-sky signal corresponding to the epoch of reioniza- 
tion, a NS orientation is preferred since this orientation has a lower antenna 
temperature over a wider LST range, and hence the sensitivity to any spectral 
signatures would be greater. 

We assume that the spectra of the sky brightness towards any direction 
are of simple power law form. At any instant, the measured sky spectrum is a 
linear combination of the power-law form sky spectra towards directions within 
the beam. This measured spectrum will no longer be a single power law and 
will require a more complex fitting form. We have computed synthetic spectra 
that would be observed with a short dipole antenna over LST and attempted 
to fit the distribution of log(T) versus log(frequency) using polynomials. The 
averaging over power law spectra results in a spectrum that is, unsurprisingly, 
dominated at the highest frequencies by the sky regions within the beam that 
have relatively flatter spectra, and is dominated at the lowest frequencies by 
sky regions that have the steepest spectral indices. Therefore, relative to the 
best fit power law model, the residuals are concave upwards in shape. We have 
attempted a fit of third order polynomials to the synthetic spectra; our fit 
takes the form: 

l°Sio T = «o + ai(logiof) + a 2 (log 10 f) 2 + a 3 (log 1Q f) 3 (12) 

where T is the temperature and / is the frequency. In order that the fitting 
function be smooth and without inflection points, the coefficients a 2 and et 3 
are constrained to be positive. The residuals of this fit to spectra that were 
synthesized at intervals of 2 hr in LST are shown in Fig. [5j and are within 
a few mK over most of the frequency range, suggesting that sky models of 
this form should suffice. It may be noted here that Pritchard and Loeb (2010) 
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Fig. 5 Residuals to fits to synthetic sky spectra. Third order polynomials of the form given 
in equation 9 were fitted to synthetic spectra. 

arrived at similar conclusions based on their fits to synthetic sky spectra that 
were derived using the global sky model of de Oliveira- Costa et al. (2008). 

The difference between sky spectra measured towards different directions 
is again essentially a linear combination of power-laws corresponding to the 
Galactic emission, since it is expected that the uniform extragalactic com- 
ponent including any re-ionization signatures would be canceled in the sky 
difference. The difference of spectra obtained during LST when the antenna 
temperature is relatively high and that when the antenna temperature is rela- 
tively low serves as a smooth calibration signal. We have divided the synthetic 
spectra computed over LST into two halves based on the median antenna tem- 
perature, computed the difference spectrum, and examined the goodness of fit 
when fitting to the form described in Equation 12. The residuals to this fit are 
within ±1 mK over most of the band. 

To summarize the discussion in this section, it is sufficient and necessary 
that the sky temperature spectrum is modeled with four parameters and we 
may adopt the form given in Equation 12. The coefficients and 03 may be 
constrained to be positive, so that the model is smooth and does not contain 
inflections. 



4 SARAS radiometer: higher order corrections for non-ideal 
behavior 

The measurement equations and calibration method described above assumes 
ideal signal propagation paths that are impedance matched at all intercon- 
nections and that there are no reflections at the interfaces. The antenna has 
also been assumed to be impedance matched to free space and to the trans- 
mission line. Practical antennas, devices, as well as transmission lines have 
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impedances that deviate from the nominal design value and this may vary 
with frequency resulting in frequency dependent reflection of signals at inter- 
faces. Additionally, when long transmission lines are involved, there may be 
impedance variations along the path owing to manufacturing tolerances and 
this may result in distributed reflections of the propagating signals. In this sec- 
tion we discuss the modeling of departures from the ideal behavior analyzed 
earlier, which arise from impedance mismatches. 

The impedance mismatch at the antenna is characterized by its complex 
reflection coefficients — also referred to as 'return loss' — which is usually fre- 
quency dependent. Hence only a part of the electromagnetic radiation incident 
on the antenna is coupled to the transmission line. This signal loss is frequency 
dependent and is over and above the direction dependent response defined by 
the normalized radiation pattern as well as any ohmic loss, which may be 
independent of the direction of incidence. 

In a real system, the LNA input impedance is complex and differs from the 
impedance of the components and transmission lines which interconnect the 
antenna and reference load to the amplifier inputs. The impedance mismatch 
at the LNA inputs are also characterized by complex reflection coefficients 
which may be functions of frequency. 

The dominant sources of noise power in the signal paths within the SARAS 
radiometer configuration are (a) the antenna temperature, (b) noise power 
from the reference load, (c) noise from the first stage of low noise amplifiers 
and, (d) when the calibration noise is on, the additional noise from that source. 
Amplifier noise has a component that adds to the input signal and propagates 
downstream along with the amplified input signal; it has a second component- 
correlated with the first — that propagates upstream against the nominal signal 
flow direction. All noise in the signal path that propagates upstream towards 
the antenna would be ideally expected to exit the system via the antenna; 
however, these are partly reflected by the mismatch at the antenna, which is 
characterized by the antenna reflection coefficient or return loss. 

In the SARAS configuration, the electronics at the antenna base are all 
packaged with modules directly connected successively or with absolute mini- 
mum lead lengths. Nevertheless, the antenna temperature, reverse propagating 
receiver noise from the two low noise amplifiers and noise power from the ref- 
erence load (and calibration noise) would undergo multiple reflections between 
the antenna, reference load and amplifier input terminals. Reflections in this 
part of the radiometer between terminations on cither side of the cross-over 
switch and power splitter give rise to spurious responses that are not canceled 
in the switching; the reflections at this stage couples power across the two arms 
of the correlation spectrometer causing the radiometer to respond to receiver 
noise as well. Multi-path propagation leads to spurious complex responses that 
appear in the real and imaginary components of the measurements as ripples 
whose period in frequency space depends on the excess traversal delay in the 
multi-path propagation. Since the spurious correlations arising from internal 
reflections across the cross-over switch and power splitter appear as complex 
responses, where as the ideal system response is purely real, the SARAS system 
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Fig. 6 Reflection coefficients of the two LNAs (Magnitude) 
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Fig. 7 Reflection coefficients of the two LNAs (Phase) 



configuration admits a modeling of the spurious additives to the real spectrum 
by modeling the ripples in the imaginary spectrum. 

We denote by symbols r a , F re f, I\ and r 2 the complex reflection coeffi- 
cients at the antenna, reference and at the inputs of the two receivers respec- 
tively. In Figs. [6] & [7] we show the measured complex reflection coefficients 
at the inputs of the SARAS low noise amplifiers in the form of plots of the 
amplitude of the reflected power, and phase of the reflected voltage, as func- 
tions of frequency. The amplifiers we have used are not a matched pair and, 
therefore, the reflectivities differ. Additionally, these functions appear smooth 
over frequency because the amplifiers are wide-band and have operating fre- 
quency range much wider than the SARAS system. In Figs. [8] & [9] we show 
corresponding plots for the measured reflection coefficient at the antenna. 

The measured complex reflection coefficients of the amplifiers and antenna 
include multi-path propagation between the measuring instrument and these 
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devices because the measuring instrument will not be a perfect match to the 
interconnecting cable. Therefore, we have used a relatively long cable for the 
measurement so that the frequency structure in the measurement that arises 
from multi-path propagation would be of relatively finer scale compared to the 
frequency structure in the reflection coefficient, and hence distinguishable in 
the measurement. The measurements have been fitted using a model for the 
frequency dependence of the reflection coefficients, which does not include — 
and hence filters out higher order terms corresponding to multiple reflections 
between the measuring instrument and device under test. It is worthwhile 
stressing that while jointly modeling the observed spectra for sky features and 
system parameters it is important to use a low-order functional form for the 
reflection coefficients and not the measured coefficients, because the measured 
coefficients are contaminated by reflections during the measurements. The fit- 
ted forms are also plotted in Figs. [6] and [7] along with measurements; deviations 
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in the form of residual ripples exist and are too small to be obvious in the plots 
shown. 

In the SARAS system the balun at the antenna is directly connected to 
the directional coupler that injects the calibration noise into the signal path; 
we assume that the dominant path length in the problem is the interconnect 
between this directional coupler and the cross over switch, which is integrated 
with the subsequent components of the receiver. Reflections over this path 
length I correspond to a round-trip phase <P — (AttIv/ f„c), where v is fre- 
quency, c is the speed of light in vacuum and /„ is the velocity factor for the 
propagation of EM signals of frequency v in the transmission line. 

Hereinafter we assume that the gain loss in signal propagation through 
the paths in the cross over switch and power splitter are equal: g S \ = gs2 = 
9Ai = 9A2 = 9= (1/v^O). 

As discussed below, we expect that there be unwanted response arising from 
multi-path propagation within the receiver between the antenna, reference and 
LNA input ports. This is predominantly because of mismatch at the antenna 
terminals; nevertheless, unwanted response also arises from reflections at the 
reference port. The reference load is a noise source followed by a precision 
23 dB attenuator, which is connected directly to the J4 port of the cross over 
switch. The impedance match at this reference port is measured to have a 
power reflection coefficient less than —60 dB. Therefore, we may expect that 
the unwanted response arising from reflections at the reference port would be 
reduced by a factor smaller than —25 dB or 0.3% of the response arising from 
reflections at the antenna port. Additionally, because the interconnect to the 
reference termination is direct and because the precision attenuators and noise 
source are wide-band devices, the response would be or low order. Based on 
the estimates made below on the magnitude of the unwanted response arising 
from reflections at the antenna port, we estimate that the reflections at the 
reference port would give rise to a smooth response that is of order 10 mK. 
We therefore neglect reflections at the reference port in the analysis presented 
here. 



4.1 Measurement equations considering reflections to first order 

In this sub-section, we derive the measurement equations for the cross pow- 
ers generated by the antenna, reference and the receiver noise voltages. We 
consider here a single reflection from the LNA inputs and the antenna termi- 
nal. Second and higher order reflections would generate higher harmonics with 
relatively smaller amplitudes and this is discussed in a later sub-section. 

We use symbols v a , v re f, v n \ and v n2 to denote, respectively, voltages 
representing the antenna noise, reference noise and the noise from the low noise 
amplifiers in arms 1 and 2 of the correlation receiver. G\ and G2 represent the 
complex voltage gains of the two receiver chains from the inputs of the low 
noise amplifiers onwards. 
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4-1.1 SARAS response to noise power from the antenna 

The calibration noise is injected into the system immediately following the 
balun via a directional coupler. Therefore, the response of the system is related 
to antenna noise and calibration noise through the same transfer function. In 
what follows we compute the response to the antenna temperature T a . When 
the calibration noise is on, T a may be replaced by T a + T cal to obtain the 
response. 

We write first the signals when the cross over switch position is '0'. The 
direct propagation of antenna signal to the low noise amplifiers results in 
voltages gv a and — gv a at their input terminals. The negative sign in the second 
term explicitly conveys the inversion in signals propagating from A to P 2 ports 
of the power splitter relative to the other paths. 

After reflection from the two receiver inputs, and from the antenna, the 
additional voltages incident at the two receiver inputs are, 

van = (A + r 2 )r a g 3 v a e 1 ' 1 ', (13) 

and 

var2 = -(ri + r 2 )r a g 3 v a e l *. (14) 

Adding the direct signals to those after reflection, the net signal voltages at 
the receiver inputs are 

Vy = g[l + (A + r 2 )r a g 2 e l4, }v a (15) 

and 

V 2 = -g[± + (A + A)AsV> a . (16) 

The response of the correlation spectrometer to these voltages is 

Aa = (G1V1) X (G* 2 V 2 *) + Pcor 

= -9 2 GiG^|[l + (A + r 2 )r a g 2 e**]\ 2 v 2 a + P cor . (17) 

At switch position T', the antenna signal travels to the receiver inputs via 
the S port of the power splitter; the reflections also travel via the S port. 
There are no signal inversions and, therefore, the response of the correlation 
spectrometer to the antenna signal in this case is 

Aa = g^GlWi + (A + r 2 )r a9 2 e^}\ 2 v 2 a + P cor . (18) 

As discussed earlier in Section 3.2, we difference the measurements recorded 
in the two switch positions to cancel the spurious correlations (P cor ) owing to 
couplings between the arms of the correlation spectrometer; this yields 

Pa = Aa — Aa 

= 2g 2 G 1 G* 2 \[l + (A + r 2 )r a g 2 e^]\ 2 v 2 a 

= 2g 2 G 1 G* 2 v 2 a [l + |(A + A)A|V + 2. 9 2 Rc{(A + AjA^}] 

= 2g 2 G 1 G* 2 v 2 a C 1 . (19) 
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In the last line above we have introduced a new symbol C\ to denote 

c y = i + 1 (a + r 2 )r„| V + 2<? 2 Re{(A + r 2 )r a e i0 }. (20) 

When calibrated for the ideal complex gain 2g 2 G\G* 2 of the correlation 
receiver, the response is wholly real and the imaginary part is expected to be 
zero. This response to antenna signal may be viewed as composed of the three 
terms, all of which are real, in the above expansion for C\. 

The first term (unity) represents the ideal case where no reflections occur. 

The second term is expected to be a slow and smooth function of fre- 
quency, essentially dependent on the reflection coefficients at the amplifier 
inputs and at the antenna. For the measured coefficients, we show the vari- 



ation of [1 + \{T\ + r2)r a \ 2 g 4 ] in Fig. 10 This second term represents the 
additional response arising from the product of signals that arrive at each of 
the receiver inputs after one reflection and is expected to be 0.05-0.4% of the 
antenna temperature or about 0.15-1.2 K. 

The third term, which is of larger magnitude than the second, depends 
on the path length I and would be a sinusoidal ripple across the observing 
band that is modulated by a relatively slowly varying function of frequency: 
the modulation amplitude is determined by 2g 2 (r\ + I2) r a . This third term 
arises as a consequence of the product of the direct propagation of the an- 
tenna signal to each of the inputs of the low noise amplifiers with the reflected 
component that arrives at the other input. In contrast to the first two terms, 
this third term is the product of signals that have a relative delay correspond- 
ing to twice the length Z, and hence has a sinusoidal variation of the form 



sva(Aiilv / f v c) whose period depends on I. In Fig. 11 we plot the net response 
|[1 + [Tx + r 2 )r a g 2 e' 1 ®} | 2 assuming that I = 5 m. The net response is domi- 
nated in its frequency structure by the ripple corresponding to the length of 
the interconnect and is expected to have an amplitude 5-10% of the antenna 
temperature, which may be as much as 25 K. The analysis illustrates why it is 
essential to have low reflection coefficients and small transmission line length 
between the antenna and receiver modules, and devise methods to calibrate 
for the spectral distortion arising from the inevitable reflections. 

It may be noted here that even if the two low noise amplifiers were matched, 
these unwanted terms would not cancel. Additionally, since all of the terms 
multiply v 2 , the response to calibration noise power would be simply the same 
expression excepting that v 2 is replaced by v 2 al . Therefore, calibration by 
dividing the antenna response by the response to calibration noise would be 
expected to remove the spurious spectral shapes in the response to antenna 
temperature, which arise from reflections. 



4-1-2 SARAS response to noise power from the reference 

The response to noise power from the reference may be computed in the same 
way as above. In this case the delayed signals arise from reflections of the 
reference noise from the inputs of the low noise amplifiers followed by reflec- 
tions from the antenna. The transfer function representing the response to the 
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Fig. 10 SARAS response to antenna temperature: the plot shows the frequency behavior 
of the second term [1 + |(i~i + -T2)-Ta| 2 9 4 ] in equation (26) that arises owing to products 
between singly reflected voltages along the two receiver chains. 
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Fig. 11 Net response to antenna temperature: the plot shows the net response 
|[1 + + r2)r a g 2 e 1 ' p ]\ 2 that accounts for signals with path delays of at most 21. We have 
assumed here that the path I corresponding to the cable interconnect is 5 m and cable 
velocity factor is 0.6 for the illustration. 



reference would be the same when the reference noise source is on as well as 
off. 

When the switch position is '0', the response of the correlation spectrometer 
to reference noise voltage v re f is 

Poref = 9 2 G 1 Gl [1 + (A - A) AgV*] [1 - (A - A) r a g 2 e^]*v 2 ref + P cor (21) 

and when the switch position is '1', 

Are/ = -g 2 G 1 G;[l + (A-A)A^][l-(A-A)A.g 2 e^]^ r 2 e/ + A ,.. (22) 
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Differencing the reference power measured in the two switch positions, 

Pref * Ire/ Ore/ 

= -2 9 2 G 1 G* 2 v 2 ref [i ~ 5 4 |(A - r 2 )r a \ 2 + i2 9 2 im{(r 1 - r 2 )r a e^}} 

= -2 3 2 G 1 G^ r 2 e/ C 2 . (23) 

As in the previous subsection, we have introduced a new symbol Ci above to 
denote 

c 2 = i - ff 4 |(A - r 2 )r a \ 2 + z2. 9 2 im{(r 1 - r 2 )r a e i0 }. (24) 

As in the case of the response to antenna temperature, the reference signal 
components arriving at each of the inputs to the amplifiers consist of directly 
propagating components — which we refer to hereinafter as d\ and d 2 — and 
components that arrive with propagation delays 21, which we hereinafter refer 
to by the symbols r\ and r 2 . The first term (which is unity) in the above ex- 
pression (Equation 24) represents the product d\d* 2 , the second term represents 
r\r 2 and the third term d\r 2 + r\d 2 . 

We assume that the response is calibrated for the ideal complex gain 
2g 2 G\G 2 . Following such a calibration, the response of the system to ref- 
erence noise power is complex unlike the response to antenna temperature, 
which was real. The first two terms in the response in equation 24 constitute 
the real part where as the third term represents the imaginary part of the 
complex response. The first term is what is expected if reflection coefficients 
were identically zero. The second term g 4 1 (A - r 2 )r a | 2 is a smooth function of 
frequency whose spectral structure depends on the reflection coefficients at the 
amplifiers and antenna. We plot in Fig.[T2]the first two terms, which constitute 
the real component of the response, for the measured coefficients in the SARAS 
system. The third term, which is imaginary, has a sinusoidal variation with 
period that depends on the path length 21 and has a modulation that depends 
on 2g 2 (r 1 - r 2 )r a . We plot this imaginary response 2 5 2 Im{(A - r 2 )r a e l<l> } 
in Fig. |13[ using the measured parameters of the SARAS receiver. 

We may make a few remarks here: 

First, in contrast to the case of the response to antenna noise, the unwanted 
response to reference noise may be eliminated by using a matched pair of low 
noise amplifiers. If J\ — T 2 , the second and third terms in Equation 24 above 
vanish. As discussed above, the corresponding terms in the response to antenna 
temperature, Equation 20, survive. 

As in the case of the response to antenna temperature, the spurious spectral 
structure arising from the d\r 2 + r\d* 2 term dominates the response. For the 
parameters of the SARAS receiver, the r\r 2 response to reference noise is 
less than 0.1% of the reference noise temperature where as the ripple in the 
spectrum owing to the d\T 2 + r\d 2 response has an amplitude that is 2-3% 
of the reference temperature. If the reference is a 300 K ambient temperature 
matched termination, the spurious ripple amplitude is as much as 9 K in 
amplitude. 

The calibration discussed in Section 3.3 envisages a measurement of the 
response to the calibration step (P ca i — -Po//), and a subsequent use of this 
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Fig. 12 Response to noise temperature of the reference: we plot here the term [1 — <7 4 |(i~i — 
r2)r a \ 2 ] that arises from products between reflected voltages along the two receiver chains. 
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Fig. 13 Response to noise temperature of the reference: we plot here the term 2<7 2 Im{(ri — 
J~2)J~ae 1 *} in the net response that represents the imaginary part. 



'bandpass calibration' as a divisor for the response of the system P Q ff so 
that the quotient would have all ripples associated with reflections of antenna 
temperature calibrated out. However, the response P ff is a difference mea- 
surement that includes the response to the antenna temperature as well as 
to the reference temperature. Even if we use matched amplifiers so that the 
response to the reference would not have any ripple across the bandpass, the 
calibration process will imprint the ripple, which would be present in the re- 
sponse to the calibration step, on to the additive response to the reference. 
Therefore, as discussed below, the analysis of the measurements of sky spectra 
must necessarily be a joint modeling of system parameters along with the an- 
tenna power, and is not a simple pipeline calibration of the measured spectra 
for additive and multiplicative responses to derive the wanted sky spectrum. 
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The differential response to antenna temperature and reference will have 
structure in the real and imaginary components dominated respectively by 
modulations with amplitudes proportional to (7\ + r 2 )r a and (-T\ — r^)/^; 
both modulations have the same period in frequency space. The differential 
response to antenna and reference temperatures may be displayed as a phase 
diagram in the complex plane as shown in Fig. 
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Fig. 14 Phase diagrams of the differential response to antenna temperature and reference 
noise. The reference is assumed to be a 300 K termination. The figure in the upper panel was 
obtained by assuming that the antenna temperature is 1000 K; the panel below corresponds 
to a model antenna temperature that is [150+200 X (//131.25) - 2 ' 55 ] K, where / is frequency 
in MHz. 



4-. 1.3 SARAS response to receiver noise 



Amplifier noise has a component that adds to the input signal and propagates 
downstream along with the amplified input signal; it has a second component — 
correlated with the first — that propagates upstream against the nominal signal 
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flow direction. A part of these upward traveling receiver noise components 
from both low noise amplifiers are reflected at the antenna (as well as at the 
reference if an impedance mismatch is encountered at that port) and would 
propagate downstream to both amplifiers and add, with a delay, to the receiver 
noise that originally propagated downstream. This causes mutual coupling of 
receiver noise between the two arms of the correlation spectrometer. 

Let /i and fa be the fractions of the noise voltages of amplifiers 1 and 
2 that propagate in reverse direction. We denote the amplifier noise voltages 
by v n i and v n2 ; the reverse propagating noise voltages are fav n \ and fav n2 . 
Herein we express the reverse propagating noise powers as fractions of the 
usual receiver noise power and the reference point is the inputs of the low 
noise amplifiers. 

In this section, we consider the usual noise powers of the amplifiers plus 
additional components arising from a single dominant reflection of backward 
propagating receiver noise, the reflection is assumed to happen at the antenna. 
At switch position '0', the net amplifier noise voltages at the two amplifier 
inputs are: 

v n ri = v nl + {fav nl - fav n2 )r a g 2 e l * , (25) 

and 

Vnr2 = V n2 + (fav n2 - fl V nl ^^V* . (26) 

The noise voltages v n \ and v n2 are uncorrelated and the correlation spectrom- 
eter gives zero mean response to terms that involve products between these 
voltages. Omitting such terms, the response of the correlation spectrometer to 
amplifier noise may be written as 



°on = (Giv nrl ) x {G* 2 v* nr2 ) + P cor 

4 n 2 r j-2„.2 , j-2 2 



= -G 1 G;[.g 4 ^{/Xi+/2X 2 } 

VK1/1* O"** + v 2 n2 far a e^}\ + P cor . (27) 

Similarly, we may write the response in switch position '1' as 

Pin = G 1 G* 2 [g 4 r 2 {f 2 v 2 nl + f 2 vl 2 } 

+g 2 {v 2 nl f*r:e->* + v 2 n2 far a e^}} + P cor . (28) 

Differencing the responses to receiver noise power in the two switch positions, 

= 2g 2 G 1 G* 2 [g 2 r 2 {f 1 vl 1 + f 2 v 2 n2 } + K^^e"** + v 2 n2 faP a e^}\ 

= 2g 2 G 1 G* 2 [g 2 r 2 .f 2 + /r^e"^]^ + 2g 2 G 1 G* 2 [g 2 r 2 f 2 + f* 2 r* a e^]v 2 n2 

= 2g 2 G 1 G* 2 [C nl v 2 nl + C n2 v 2 n2 }. (29) 

Once again we introduce new symbols C n i and C n2 that are 

c nl = [(ghr a f + [hr a e l *y] (30) 

and 

C n2 = l( 9 f 2 r a ) 2 + (farae 1 *)]. (31) 
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Following the notation adopted in previous subsections, the first terms in 
Equations 30 & 31 are r\r 2 type terms and the later are d\r 2 and T\d 2 type 
terms. The first terms are real and subdominant; they are smoothly varying 
functions of frequency that depend on the frequency dependance of the receiver 
noise, that of fractions /i and / 2 , and that of r a . 

Assuming that the fractions /i and / 2 are real quantities, and writing 
r a = |-Ta|e~* a , the second terms in Equations 30 & 31 may be expanded 
together as follows: 

= vlj 1 \r a \e-^+^ + vlj 2 \r a \e^+^ 
= + M 2 l2 )cos(<z> + <P a ) - 

i(fiv 2 nl -f 2 v 2 n2 )sm^ + ^ a )}. (32) 

This is the dominant term representing the summation of d\r 2 and r\d 2 type 
terms and this part of the response is complex. The real and imaginary parts 
are both modulated sinusoids with tt/2 phase difference, whose amplitudes are 
proportional to (fiv^ + f2V n2 ) an d (h v ni ~ hv n2 )- I n the complex plane, the 
trace of the response is ellipses with major and minor axes l^aK/i^i + f2V 2 l2 ) 
and lAK/iw^ - / 2 v 2 2 ) respectively. 

We have experimentally measured the fraction of the receiver noise that 
travels upwards in the SARAS system. To reconfigure for this measurement, 
we first remove the antenna as well as the switch and power-splitter assembly. 
While one LNA input is terminated with a matched impedance, the LNA 
under test has a cable of length I = 5 m connecting its input directly to 
the output of the directional coupler. The main line input of the directional 
coupler is terminated in an 'open' load and the coupled port is terminated with 
a matched impedance. The auto-correlation spectrum of the output of the LNA 
under test is recorded. In this configuration the fractional noise from the LNA, 
which travels upstream, suffers total internal reflection from the open load. 
Therefore, the measured auto-correlation has contributions from the receiver 
noise that travels directly downstream, and a contribution arising from the 
additional correlations owing to the reflected receiver noise. Since the direct 
and the reflected receiver noise voltages have a phase difference corresponding 
to the path difference of 21, the correlation spectrum is oscillatory in nature. 
We write below the analysis for the LNA with label '1'. Following the same 
notation used above the net noise voltage at the input of the LNA, considering 
one reflection of the upward traveling receiver noise at the open load, is given 
by: 

Vnrl = V n i + flV nl e L ' t ' . (33) 

The power spectrum corresponding to the noise voltage at the LNA may be 
written as: 

Pnl = {Giv nrl ) x (G\v* nrl ) 

= G^ ? 2 ll (l + / 1 2 + 2Re{/ 1 e^}). (34) 
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The first term in parenthesis (unity) represents the autocorrelation of the 
receiver noise traveling directly downstream: d\d\. The second term is the 
autocorrelation of the upward traveling receiver noise that is reflected back 
from the open load: r\r\. The third term represents the correlation {d\r\ + 
d\ri) between the downward traveling receiver noise and its upward traveling 
component, following reflection. 

Bandpass calibration of the above measurement is done in a modified setup. 
The 5-m cable is omitted, the 'open' load is replaced with a matched termi- 
nation, and a calibration noise source is connected to the coupled port. The 
calibration noise is toggled between 'on' and 'off' every 0.7 s and the band-pass 
calibration is derived by subtracting the recorded power spectra in 'off' state 
from that in 'on' state. The band-pass calibration is given by: 



cal 



= G\vt al . (35) 



Dividing P n \ by P ca i and multiplying by the temperature T n \ of the calibration 
noise source, the calibrated spectrum of the power from receiver 1 is given by: 

Tn = T„i(l + H + 2Re(/ie li) )). (36) 

Higher order reflections of LNA amplifier noise between the LNA input 
and the open end of the transmission line would add signals with delays that 
are integral multiples of 21, These would generate the ripples in the measured 
spectrum with higher harmonics, amplitudes of which are successively reduced 
by about 10%. We estimate the fractions fx and f 2 via fits to the measured 
spectra for only the fundamental; therefore, we have neglected writing the 
analysis above as a series expansion representing the higher order reflections. 

We fit the measured data to the form ao + aisin(a2/), where ao represents 
the mean receiver noise power and a\ represents the amplitude of the ripple. 
The residual to the fit contains higher order reflections as well as a residual 
ripple due to the frequency dependence of cable loss and reflection coefficient 
at the LNA input. Following calibration, the measured spectra corresponding 



to the two receivers are shown below in Fig. 15 along with the fits. We thus 
obtain /i = 9.9% and f 2 = 4.67%. 

We plot g 2 Tl{fl + /f} for the SARAS system in Fig. Il6] adopting the 



above measured values for the fractions fi and f 2 - The sum of these terms in 
Equations 30 & 31 is expected to manifest in a response with magnitude less 
than 10 mK for receiver noise temperatures of a few tens of Kelvin. For ( = 5m, 



we plot in Fig. 17 the real and imaginary components |-Ta|(/i + /2)cos(^ + ^ a ) 
and l-TaK/i — / 2 )si n (^ + ^a) in the upper panel. The complex response is also 
shown as a phase diagram in the complex plane in the lower panel of the figure. 

The response may have an amplitude as much as 5% of the receiver noise 
temperature. The imaginary response is always less than the real and is offset 
ir/2 in phase. If the amplifiers form a matched pair and /i = f% and = 
(Vrfi), the imaginary term vanishes. 
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Fig. 15 Calibrated spectra of the power from the receivers with 5 m cables connected to 
their inputs; the cables are open at their ends to reflect the amplifier noise traveling out of 
the input ports. Also shown arc fits using sinusoidal forms. The upper panel is for the first 
receiver chain and the lower panel is for the second chain. The y-axis is in arbitrary units. 



4.2 Measurement equations considering higher order reflections within the 
system 

The analysis described above yields the response of the SARAS spectrometer 
considering only first order reflections. We have therein included effects aris- 
ing from a single reflection at the antenna. However, precision measurements 
of the cosmic radio background with realistic systems, which aim to detect 
weak features from the epoch of reionization, would require a more complete 
analysis that takes into account multiple or higher order reflections of the sys- 
tem temperature components. The higher order reflections generate ripples of 
higher order in the spectra measured by the correlation spectrometer; these 
would be at harmonics of the fundamental ripple discussed in previous sec- 
tions. Herein we develop the measurement equations that include effects of 
higher order reflections in the form of series expansions. 
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Fig. 16 Response to receiver noise: the figure shows the term g 2 T^ {/j* + /f } in Equation 35 
that arises from r\T^ products between reflected signals of upward traveling receiver noise 
voltages from the two low noise amplifiers. 
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Fig. 17 The complex response to receiver noise; what is plotted is the sum of factors that 
multiply the receiver noise to yield the response. On the top panel is shown the real and 
imaginary components of the response; in the lower panel the response is shown in the 
complex plane. 
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4-2.1 SARAS response to noise power from the antenna 

We begin by considering the response to the signal from the antenna. Con- 
sidering the state of the system with cross over switch in position '0', the net 
signal voltages at the pair of amplifier inputs are 

V l = g[l + (A + r 2 )r a g 2 e^ + (A + r 2 ) 2 r 2 g 4 e i2f + 

+ (A + A;"/;;'/r 'v'"" + h (37) 

and 

v 2 = -g[i + (A + AOAsV* + (A + rrfrlg^ 2 * + 



+ (A + r 2 ) n r?g 2 * n e in * + }v a . (38) 

Using the symbol T to denote the expression (A + A)Aff 2 , 

Vi = g[l + Ye 1 * + TV 2 * + .... + r n e in * + }v a 

oo 

= . 9 (]Tr"V m > a (39) 

and, similarly, 

oo 

V 2 = -g{Y J rm z im *W (4°) 

m=0 

The net response of the correlation spectrometer for this switch position may 
be written as 

P 0a = (GlVl) x (G 2 V 2 )*+P cor 

OO oo 

= -g^G^T^^T^^yvl + P cor 

oo 

= -g 2 GiG;[^ \r m \ 2 (l + Re(Te^) + Re(T 2 e 



,2<1>\ 



m— 



+.... + Rc(T"e mS> ) + )K+P cor . (41) 

At switch position '1', {P\ a — P cor ) — —{Poa — Pcor)- Subtracting the cross- 
power in the two switch positions to cancel P cor , we get 

Pa = P\a ~ Poa 

oo oo 

= 2.g 2 G 1 G^[]r (|T m | 2 ) ]T Re(T k e^)]vl 

oo 



2g 2 G 1 G* 2 [J2(\T m \ 2 )]vl + 

oo 

2 3 2 G 1 G^[Rc(re^)^(|r m | 2 )]« 2 + 



m=0 
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2g 2 G 1 G* 2 {Rc(r 2 e* 2 *) ^(|T"f )}v 2 a + 

oo 

+ 2g 2 G 1 G* 2 [Rc(T n e m4 ') ^ (\T m \ 2 )]v 2 a + (42) 

m=0 



The terms in the above expansion may be interpreted as follows. The cross 
correlation of signals arriving directly at the two LNA inputs, as well as of 
signals arriving at the pair of inputs after having undergone identical number 
of reflections, would result in spectral responses that would vary smoothly as 
a function of frequency: these are the terms in the first line of the expansion. 
However, every additional path delay of 21 in which the signals reflect at an 
LNA input and antenna cause the incident voltage amplitude at the LNA 
input to reduce by factor T = (Ji + r 2 )r a g 2 , which is at most 5%. The cross 
correlation response to components that have undergone successive reflections 
reduces by factor T 2 = ((r\ + r 2 )r a g 2 ) 2 , which is about 0.3% of the response 
to the adjacent lower order reflection. Beyond the first order reflection, high 
order terms yield responses below about 1 mK; therefore, we may neglect their 
contributions to the net response of the system. 

The second line in the expansion represents products between antenna 
signals that are reflected p and p + 1 times. The correlation between any 
two successive reflected components would result in a spectral ripple with the 
fundamental period. However, in the response the amplitudes of such ripples 
in successive terms would again be reduced by factors Y 2 . The first term in 
this series may be as large as 25 K, the second term would be at most 75 mK 
in amplitude and later terms less than 1 mK. 

The third line in the expansion represents the first harmonic, which arises 
from correlations between antenna signals arriving at the LNA inputs with 
relative delays of 4L The first term with this first harmonic ripple period 
merits consideration since its amplitude is only a factor \Y\ (which could be 
as much as 5%) below the amplitude of the fundamental. The second term is 
an additional factor T 2 lower and hence would not be more than a few mK. 

The correlation between direct component along one path with the re- 
flected component along the other path which has undergone 'rC reflections 
would result in a ripple which is (n — l)'th harmonic of the fundamental ripple 
frequency. This (n— l)'th harmonic has an amplitude that is a factor |7 , | n_1 of 
the fundamental. In the SARAS system, the third and higher harmonics would 
have amplitudes below a few mK and may be neglected while modeling EoR 
measurements. Thus, for the modeling of SARAS spectral measurements that 
are aimed at detecting the all-sky EoR signatures, the response to antenna 
noise power may be approximated as: 

P a « 2g 2 G 1 G* 2 [l + \T\ 2 +Rc{Te^)(l + \T\ 2 ) + Rc(rV 2<f )(1 + \r\ 2 )]v 2 a . (43) 
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4-2.2 SARAS response to noise power from the reference 

In this subsection we analyze the response of the SARAS system to noise 
temperature of the reference. As above, we consider only reflections from the 
antenna terminal and neglect reflections from the reference terminals where the 
impedance of the termination is adequately matched. In case of the antenna 
temperature, successive reflections of the signal from the antenna to the LNAs 
traverse the same path as that of the direct incidence. However, in the case of 
the reference noise voltage the multiply reflected signals travel paths between 
the antenna and LNA inputs and not the path of direct incidence between 
the reference and LNA inputs. This asymmetry results in that the response 
to reference noise power is dominated by a Fourier series of ripples in the 
imaginary component of the complex response, unlike the case for the response 
to antenna noise power that is dominated by a Fourier series of ripples in the 
real part. 

We begin by considering the response of the system with cross over switch 
in position '0'. The notations used below are the same as in Section 4.2.1. The 
net signal voltages at the pair of amplifier inputs are 

V 1 = g[l + (A - r 2 )r a g 2 {e^ + (A + r 2 )r a g 2 e 2i4 ' + (A + A) 2 AVe l3 * + 
+ (A + A) ( ™~ 1) A^~ 1) 5 2( "~ 1) e m<l> + }]v ref (44) 

and 

V 2 = g[l - (A - A)A.g 2 {e jS ' + (A + A)A.gV 2<f + (A + A) 2 AVe l3<f + 

+ (A + J T 2 )(™-i)A<"- 1) ff 2(rl - 1) e m<1, + }]v ref . (45) 

We continue to use the symbol T to denote (A + A)A.9 2 , introduce an ad- 
ditional symbol <F to denote (A — A) A3 2 , and reduce the above equations 
to 

V 1 = g[l + ^{e^ + Te lM> + TV 3 * + .... + T {n - V > ] e m<p + 

OO 

= #+^{E r(m ~ 1)e ""^Ke/ 
rn — 1 

and, similarly, 

00 

V 2 = g[l - r ( " , - 1) e <m *}]tve / . (47) 

m— 1 

The net response of the correlation spectrometer for this switch position may 
be written as 

Are/ - (GlVi) X (G 2 V 2 )*+P cor 

OO OO 

= g 2 G x Gl[l + y(Y j r^e^)][l -W^T^e^)]*^ + P cor 

p=l q=l 



}}v re f 

(46) 
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OO OO 

= g 2 G 1 G* 2 [l- |^| 2 (^T( p - 1 )e ip *)(^r( g - 1 ) e i<z *)* 

p=l q=l 

oo 

+2ilm{#( J2 T^-^e™*)}^ + P cor . (48) 

m— 1 

As before, the responses in switch positions '0' and '1' have the relationship: 
(Piref — Pcor) = — (Poref ~ Pcor) ■ Subtracting the cross-power in the two switch 
positions to cancel P CO r, we get 

Pref — Plref Porcf 

= -2g 2 G\G 2 v 2 ef 

OO 

+2g 2 G 1 G* 2 \&\ 2 ^ l^™" 1 ! 2 ! 1 + Re(Te 1 *) + Re(rV M> ) + 

771—1 

+ Rc(T n - 1 e l< - n - 1 ^) + }v 2 ef 

OO 

+i4g 2 G 1 G* 2 lm{V(Y^ T m ~ 1 e im<p )}v 2 ef . (49) 

m—l 

Compared to the response to antenna temperature, the response to ref- 
erence noise is complex. The response to multi-path propagation of reference 
noise is less than that due to multi-path propagation of antenna noise owing 
to the appearance of an additional factor \P. This factor ]P = (A — P 2 )P a g 2 
has a magnitude of about 2% in the SARAS system. 

The first term in Equation 49 represents the nominal system response 
neglecting multi-path propagation, which is the response to reference noise 
traveling directly to the pair of LNAs in the two arms of the correlation spec- 
trometer. When calibrated, this term corresponds to the ambient temperature 
of the matched termination provided by the reference, which is about 300 K. 

The second term is a 'real' response that is similar in mathematical descrip- 
tion to the series of terms in Equation 42. However, because of the term \ip\ 2 
that appears as a common multiplier, all terms corresponding to differential 
delays exceeding 21 may be neglected since the response to these are expected 
to be less than 1 mK. Additionally, only the m — l term in the expansion 
need be considered. The 'real' response will, therefore, have a spectral ripple 
with amplitude less than 6 mK and the ripple will have the fundamental pe- 
riod corresponding to a differential delay of 21; harmonics will have amplitudes 
well below 1 mK. 

The third term in Equation 49 represents the 'imaginary' response. This is 
terms corresponding to correlation products between the direct propagation 
of reference noise along one arm with the rt'th reflected component along the 
other path. These products represent terms of type d\r 2 and r\d 2 and the 
net correlation is imaginary. The imaginary response is a Fourier series of 
ripples with zero expectation for the mean response. In the SARAS system, 
the fundamental is expected to have amplitude about 6 K; harmonics have 
their amplitudes reduced progressively by factors about 5% and, therefore, we 
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may limit consideration to just the first and second harmonics while modeling 
the spectra to 1 mK accuracy. 

Retaining the significant terms reduces Equation 49 to 

P ref w -2 ff 2 GiG;[l - \W | 2 {1 + Re(Te^)} - 2tlm(!?e < * + VTe t2 * 

+*r 2 e**))v 2 ref . (50) 

4-2.3 SARAS response to receiver noise 

Components of the receiver noise traveling from the LNA inputs opposite to 
the signal flow path encounter the antenna and reference terminations. As 
discussed earlier we neglect the reflections from the reference because of the 
excellent impedance match provided at that port and consider only the dom- 
inant reflections from the antenna terminal. In this section, we develop the 
expressions for the net complex response of the SARAS correlation spectrom- 
eter to receiver noise; this analysis considers multiple reflections of the LNA 
noise between the amplifier inputs and antenna terminals. When the cross-over 
switch connects the antenna to the E port of the power splitter the sum of 
the upward traveling noise voltages from the two LNAs reaches the antenna 
and when the antenna is connected to the A port of the power splitter the 
difference between the voltages propagates to the antenna. Following similar 
notations as in previous sections and considering the state of the system with 
cross over switch in position '0', the net signal voltages at the pair of amplifier 
inputs are: 

Vi = vm + (fivm - f 2 v n2 )r a g 2 e^ + {f lVnl - / 2 ^ 2 )(A + r 2 ) {T a g 2 f e i2 * 

+ + (hv n i - f2v n2 ){r l + r 2 )("- 1 )(r aff 2 )"e m<s + 

OO 

= vm + (/i«„i - f 2 v n2 ) (A.9 2 ) m (A + r 2 ) (m_1) e im * (51) 

m— 1 

and, similarly, 

oo 

V 2 = v n2 - (fiv nl - f 2 v n2 ) (Afl 2 ) ro (A + r 2 )( m -Ve im *. (52) 

m—1 

The complex response of the correlation spectrometer is: 

P n = (Gl^l) X (G 2 V 2 )*+P cor 

= GiG* 2 [v nl v* n2 + {fiv nl - ! 2 v n2 )v* n2 x ~ (/l«nl - f 2 Vn 2 )*V nl X* 

-(/l«nl - f 2 V n2 )(.flV nl - f 2 V n2 )*XX*} + Poor, (53) 

where we have used the symbol \ to denote the series expansion: 

oo 

x = E (A<? 2 ) m (A + r 2 )( m ~v e im * 

m—1 

oo 

= g 2 P a e li r m - 1 e*( m - 1 >*. (54) 

m—1 
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The time integral of the product of noise voltages from the two LNAs is 
expected to vanish; therefore, the SARAS response only contains the terms 

P 0n = -G1G2 [VniiflX* + fhx*) + V 2 n2 (f2X + fhx*)} + Poor- (55) 

With cross over switch in position '1', a similar analysis leads to the re- 
sponse 

P ln = G 1 G* 2 [vl 1 (f lX * + fhx*) + V 2 n2 (f2X + fhx*)} + Poor- (56) 

Differencing the measured complex responses in the two switch positions, 

Pn Pin 

= 2G 1 G* 2 [v 2 nl (f lX * + fhx*) + v 2 n2 (f 2X + fhx*)} 
= 2G 1 G* 2 [v 2 nl f lX * + v 2 n2 f 2X + (v 2 nl f 2 + v 2 n2 f 2 )xx*} 

OO 

= 2g 2 G 1 G* 2 f 1 v 2 nl r:e-^(Yl r^ 1 ^™" 1 )*)* 

m — 1 

oo 

+2 g 2 G 1 G* 2 f 2 v 2 n2 r a e^ r™" 1 ^" 1 " 

m— 1 

oo oo 

+2.g 2 G 1 G;/ 1 2 « I 2 ll |r a | 2 (^TP- 1 e l (f- 1 ) <f )(^T"- 1 e ^«- 1 )^)* 

p=l q=l 

CO oo 

+2.g 2 G 1 Gl/ 2 « 2 2 |r a | 2 (^TP- 1 e l (f- 1 ) <f )(^T«- 1 e ^«- 1 ) <f )*. (57) 

p=l g=l 

The first line in the above expansion represents correlation products be- 
tween the receiver noise (y n i) of the LNA in Arm 1 with the component of 
this receiver noise (fiv n i) that propagates in the reverse direction from the 
input port of that LNA and arrives at the correlator along Arm 2. This com- 
ponent may incur multiple reflections between the LNA inputs and antenna 
and hence the response would be a harmonic series of ripples in the spectral 
response. The second line in the above expansion represents a similar response 
for the receiver noise of the LNA in Arm 2. The third line represents corre- 
lation products between the reverse propagating receiver noise from the LNA 
in Arm 1, which arrives at the correlator along both arms after multiple re- 
flections between the antenna and LNA terminals; these signals have relative 
delays of zero and multiples of 21. The fourth line is a similar response to 
reverse propagating receiver noise from the LNA in Arm 2. 

We adopt parameters relevant to the SARAS system, where fi and f 2 are 
in the range 5-10%, and the magnitudes of Y and r a are, respectively below 
5 and 10%. The receiver noise temperatures are about 115 K. 

The responses written as the third and fourth lines in the above Equa- 
tion 57 contain terms that arise from receiver noise that propagates along the 
two arms after multiple reflections at the antenna terminals. This response 
is 'real' and with magnitude less than 0.01% of the receiver temperature, or 
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about 10 mK. Only terms corresponding to singly-reflected receiver noise prop- 
agating along the two arms are significant. These are the terms corresponding 
to p = 1 and q = 1 and for these terms the response is relatively a smooth 
function of frequency 

All ripples in the response that are above mK in amplitude arise from the 
first two lines of the expansion. The ripple with fundamental period corre- 
sponding to path delay of 21 is expected to have an amplitude of less than 
one per cent of the receiver temperature, or less than 1 K. Successive har- 
monics have amplitudes reduced by \T\, which is less than 5%. Therefore, we 
may consider the first harmonic alone while modeling observations for EoR 
signatures. 

Retaining these significant terms, the SARAS response to receiver noise 
may be approximated to: 



If the two LNAs form a matched pair then the imaginary component of the re- 
sponse would vanish. In SARAS, we expect a 3:1 ratio between the amplitudes 
of the sinusoidal response in the real and imaginary components. 



4.3 Bandpass Calibration revisited 

The primary calibration of the complex gain of the correlation spectrometer 
is based on the calibration noise source, which is switched on and off in each 
of the switch positions. This noise source is coupled to the signal path via 
a directional coupler, and the switching of the noise power is achieved by 
turning the noise source on and off rather than switching the coupling of the 
power to the signal path. Both these approaches are aimed at maintaining the 
impedance of the system interconnects unchanged during the calibration cycle, 
so that the calibration signal suffers the same reflections within the system as 
the antenna signal, and the calibration may be used to calibrate the spectral 
response to the antenna signal including effects of multi-path propagation. 

Following the same notation as in Section 3.2, we denote the responses of 
the system when the calibration noise source is in 'on' state as P\ ca i and Pocai 
in the two states of the cross over switch. The responses when the calibration 
noise source is 'off' are Pioff and Pooff- 

The net response of the system is: 



P n » 2g 2 G 1 G* 2 [(f 2 v 2 nl + f 2 v 2 n2 )\r a \ 2 + 

(fiv 2 nl + hv 2 n2 ) x Re{r o (e i * + TV 2 *)} + 
i(fiv 2 nl hv 2 n2 ) x Im{r a (e^ + Te i2 *)}]. 



(58) 



Poff = Ploff ~ Pooff 

= (Pla — Poa) + (Plref ~ Poref) + (-Pin _ Pon) 
= (Pa + Pref + Pn) 



(59) 
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The double difference: 

Pbpc — Peal Poff 

= (Plcal - Pocal) - (Ploff - Paoff) (60) 

represents the response of the system to calibration noise alone and has the 
same form as the response to antenna noise power: the response is given by 
Equation 42 excepting that the antenna signal v\ is replaced by v 2 cal . The 
response to calibration noise is wholly 'real'. The bandpass calibration is 

oo oo 

P bpc = 2. 9 2 G 1 G"[£(|rH 2 )]TRe(rV^)K 2 ai - 

oo 

= 2g"G 1 G* 2 [Y / (\T m \ 2 )y cal + 

oo 

2g 2 G 1 G* 2 [Rc(Te^)Y / (\r m \ 2 Mai + 

m—0 

oo 

2. 9 2 G 1 G;[Re(TV 2 *) £ (|TH>L + 

oo 

+ 2 3 2 G 1 G^[Re(T" e ^) ]T {\T m \ 2 ))v 2 cal + (61) 

m=0 

Bandpass calibration of the net response is complex arithmetic and in- 
volves dividing the net response of the system (Equation 59) with the above 
Pbpc so that the complex frequency dependent gain 2g 2 GiG2 is calibrated 
for. Additionally, the measured spectrum is also given a temperature scale 
by multiplying the quotient with the strength of the calibration noise T ca i in 
temperature units: 

rpl (Pa + Pref + Pn) rp rp . \Pref + Pn-\rp i co \ 

L a = p * J-cal = J-a + l 5 \J-cal- (62) 

*bpc +bpc 

The bandpass calibration yields a complex result. The real part is the sum 
of the antenna temperature spectrum and bandpass calibrated (P re f + P n )- 
The SARAS system requires a joint modeling of the real and imaginary parts 
of T' a to separate the contributions from reference noise and receiver noise 
and isolate T a . Useful in this modeling is the expectation that T a has no 
contribution to the imaginary part of T' a . 

5 Modeling of SARAS data 

In this final section, we discuss the data products from the spectrometer and 
an approach to modeling the system towards deriving useful measurements of 
the cosmic radio background and of re-ionization signatures embedded therein. 
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Fig. 18 shows a diagram of the SARAS spectrometer depicting the noise pow- 
ers and corresponding signals that are of significance: noise powers from the 
antenna (P a ), reference source (P re f) and the receiver (P n ) and corresponding 
voltages v a , v re f and v n i and v n i. Also depicted are the gains g, G\ and G2 
within the power splitter and in the two arms containing the receiver chains, 
the voltage reflection coefficients r a , Zi and I2 at the antenna terminals and 
at the two LNA inputs, and the noise voltages originating in the LNAs in the 
forward (w nl and v n i) and reverse (/iV„i and f2V n 2) directions. The figure 
also shows the equations for the responses of the system in terms of C\ , C2 , 
C„i and C' n 2, which are functions of many of the above system parameters 
depicted here. 



5.1 The SARAS measurement set 

During an observation, SARAS measures two complex bandpass calibrated 
cross power spectra. One of them is a differential measurement of the sky tem- 
perature with reference to a 'cold' reference load and the other is the difference 
between the sky temperature and a 'hot' reference. We have used the terms 
'hot' and 'cold' reference loads to refer to the on and off states of the refer- 
ence noise source. The receiver noise also contributes to the differential output 
measurement due to cross coupling between the receiver channels, which is 
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primarily a consequence of the mismatch at the antenna terminals. As dis- 
cussed previously, in these differential measurements the sky power (antenna 
temperature) appears only in the real part. Following on from Equation 62, 
the above differential measurements (T a0 and T' al respectively) may be written 
in temperature units as: 

TaO = G a T s ky — -prT re f + -^-T nl + — ^-T n2 (63) 
Ol Oi Oi 

and C C C 

T 'al = G a T sky - -^-T re fi + -j^-T nl + -^-T n2 , (64) 

where T a , T re f 0} Tre/i, T n \ and T n2 are the antenna temperature, hot and cold 
reference load temperatures and the two receiver noise temperatures respec- 
tively. G a is the antenna gain, which is a function of frequency and defines the 
antenna bandpass. The frequency response of the antenna includes the change 
in gain of the fat dipole versus frequency as well as the frequency dependence 
in the response to sky temperature that arises from the ground reflection; both 
these are multiplicative and modeled as a single antenna gain term G a . The 
primary data product or measurement set from SARAS is a series of T a0 and 
T al complex spectra in the frequency range 87.5 to 175 MHz with a spectral 
resolution of 85.4 kHz; a pair is computed every 4.2 s. The SARAS antenna 
has a frequency independent pattern that is a dipole; the axis of the dipole is 
parallel to the ground and hence half the pattern is on the sky and the other 
half sees an absorber-covered ground. 

From the measurement set additional calibration products are derived: one 
dependent only on the reference temperature step and the other on differential 
sky temperature. The first is obtained by differencing the two measurements 
made with the reference termination in the 'hot' and 'cold' states, i.e., by 
differencing Equations 63 & 64 at each integration time to get 

C 2 

T a0 - T al = ^riTrefl - T ref0 ). (65) 

To derive the second, we difference sky spectra measured at different local 
sidereal times. For example, a particular realization of this differential sky 
spectrum may be obtained by dividing the measurement set into two equal 
halves based on the median antenna temperature during the entire observation, 
the measurements with above-median antenna temperatures form the 'hot' sky 
spectra and the other half forms the 'cold' sky spectra. We may write the mean 
antenna temperatures over 'hot' and 'cold' skies as T a h = G a Thot sky and T ac = 
G a T co id sky respectively in terms of the mean sky brightness temperatures over 
the 'hot' and 'cold' skies. We compute the difference between the hot and cold 
sky spectra to get 

T a h ~ T ac = Tah — T ac = G a {Thot sky — T co id sky)- (66) 

The first calibration product is a measure of the response to reference load, 
including calibration. The second yields calibration products that serve for 
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joint modeling of the parameters for the foreground sky model and those of 
G a , the frequency response of the antenna or its bandpass shape. 

Critical to bandpass calibration of the measurement is the calibration noise 
source, which is injected into the signal path via the directional coupler at the 
antenna. The calibration of the equivalent noise temperature and flatness of 
this noise source is part of the observing strategy; this involves replacing the 
antenna with a matched load with temperature sensor and thermally cycling 
the termination and analyzing the measurement set to solve for the charac- 
teristics of the calibration noise source. Unfortunately at the precision of the 
SARAS system, the complex impedance of terminations appear to change with 
temperature, and this manifests as a scaling and shift in the measured ripples 
arising from reflections of reference noise and receiver noise at the termination. 
The process of primary calibration includes solving for this shift and scale in 
the imaginary part of the spectrum and correcting for this effect. 

5.2 Instrument model 

To the degree of complexity that we have considered in our analysis presented 
here, the SARAS system is described in terms of eight system characteristics. 
The measurement set may be related to the sky spectrum via the measurement 
equations, which depends on the gain of the antenna (G a ) and that of the 
power splitter (g), complex reflection coefficients at the antenna terminals 
(r a ) and at the inputs of the low noise amplifiers (A & r 2 ), the time delay in 
the propagation path between the antenna and receiver (r = <£/(27r/)), and 
the fractional LNA noise voltage that propagates out the LNA input terminals 
(/i & / 2 ). Of these, t, g, fi and f 2 are modeled as independent of frequency. 
In expanded form the defining model for the instrument is in Equations 42, 
49 & 57. 

As discussed in previous sections, the SARAS system is configured in ex- 
perimental setups to measure r a , 7\, I2, /1 and / 2 . Antenna measurements 
yield the absolute scaling factor G a at the band center; an estimate for r is 
obtained from the length and characteristics of the coaxial cable interconnect 
between the antenna and receiver. All of these yield good initial guesses for 
the system parameters, and fits to these measurements provide fitting forms 
whose parameters may be refined in the modeling of the data. 

5.3 Sky model 

The antenna temperature T a is the product of sky temperature T s k y and the 
antenna gain-loss function G a . Therefore modeling the antenna temperature 
involves modeling the sky spectrum and G a . 

At any instant T s k y is the average of the sky brightness temperature 
weighted by the antenna beam; the sky brightness along any line of sight 
is the cosmic radio foreground (Tf g ) due to extragalactic radio sources and 
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the Galactic synchrotron emission in addition to redshifted 21 cm background 
(Tb g ) from cosmological distribution of neutral hydrogen gas. The foregrounds 
are modeled using a functional form similar to the one used by Pritchard and 
Loeb (2010). As we have shown in Section 3.4, to reduce the modeling resid- 
uals to mK levels the beam-averaged sky spectrum requires modeling with 
four parameters and using a third order polynomial in log-tcmpcrature versus 
log-frequency space. The foreground sky spectrum is described using the form 

log 10 T /9 = a + Oi(log 10 (-f )) + a 2 (log 10 (^-)) 2 + a 3 (log 10 (f )) 3 , (67) 

Jo Jo Jo 

where / c =131.25 MHz is the mid frequency of our band. 

The simplest model for the background reionization feature that may be 
expected imprinted on the cosmic radio background is a toy model 'step'. 
This has a tanh-form discontinuity that is parameterized by the redshift of 
reionization z r and the duration of reionization A z (this model was adopted 
in the WMAP7 analysis; Larson et al. (2010)): 

T bg {z) = ^(/^It-M^f 1 ) + !]• (68) 

T21 is the amplitude of the reionization signal. More realistic models for reion- 
ization lead to expectations of the forms derived in Pritchard and Loeb (2008) 
and in Pritchard and Loeb (2010), whose detection requires band widths in 
excess of an octave and ability to discern broad spectral signatures in the 
presence of orders of magnitude more intense continuum foregrounds. 



5.4 Modeling the primary and the derived spectra: a hierarchical approach 

The primary measurement is Equation |63[ in which the antenna temperature 
appears in the real component. The imaginary part of Equation |63| and the 
derived calibration products described in Equations [65] & [66] do not contain 
features of the uniform cosmological radio background and reionization signa- 
tures: the modeling of these spectra calibrate the instrument parameters. 

The first step in the modeling of the measurement set is the modeling of 
the difference spectrum computed between 'hot' and 'cold' reference, which is 



described in Equation 65 This depends on parameters T = (I\ + r 2 )r a g 2 and 
<P = (A — r2)r a g 2 , whose complex reflection coefficients are measured in the 
laboratory except for the gain term g. This also depends on the propagation 
delay r and fitting for t is expedited by having a long cable length for the 
interconnect between the antenna and receiver, which would cause a number 
of cycles of ripple in the spectrum. Such an approach is preferred when low 
order background spectral features are sought to be detected as, for example, 
in measurements of the absolute spectrum of the cosmic radio background. 
However, in observations aiming to model the data for plausible 'steps' arising 
from cosmological reionization, a short interconnect is preferred. 
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This step of the modeling, which has a focus on T and If - , underscores a 
positive aspect of the SARAS design. The modeling at this stage is of the ratio 
of parameters C2 and C\. These parameters include ripples across frequency 
space of form Ye 1 ® and ^e**, which are modulated by the very terms T and 
H? that are also parts of additive functions in Ci and C\ that manifest as low 
order spectral features. The design has the advantage of allowing solutions for 
the low order terms via solving for the modulations in the high order ripples 
in the spectrum, which may be introduced by choice of cable length to be of 
a high order and disjoint in Fourier space from astrophysical features. 

The next step in the hierarchical modeling is to expand the modeling to 
include the imaginary component of the primary measurement Equation |63| 
This step introduces instrument parameters related to receiver noise, this sec- 
ond step adds the term (/iv, 2 ii — H v ^i2l to the instrument model. 

The last step in this hierarchical modeling is to include the real component 
of the primary measurement together with the calibration product that is the 
difference between 'hot' and 'cold' sky: Equation [66] Additional instrument 
parameters related to receiver noise are now included in the model: {fiv^i + 
h v 1i2) an d (/1 w ni+/l v n2)- The relative gain of the antenna versus frequency — 
G a — is another parameter that is introduced into the solver at this stage and 



is constrained by the difference Equation 66 



The SARAS system design thus allows a hierarchical approach where the 
terms defining the instrument model may be solved for first in the imaginary 
component of the measurement set and in calibration products, which do not 
have any response to spectral features in the cosmic radio background, before 
moving on to the final stages of modeling wherein terms defining the sky model 
are brought into the solver. Modeling at all stages may follow the Metropolis- 
Hastings algorithm in which all the unknown parameters are constrained si- 
multaneously within the Bayesian framework of analysis. The approach also 
allows for the parameters describing the reflection coefficients r a , Ji and r 2 
to be refined in the solver. A detailed discussion on the Monte Carlo Markov 
chain (MCMC) methodology for estimation of cosmological parameters and 
simultaneous solution for instrument coefficients can be found in Harker et. 
al. (2012). 



6 Postlude 

Analysis of SARAS measurements for deriving constraints on parameters defin- 
ing the sky spectrum and its features requires an analysis pipeline that derives 
likelihood estimates for sky model parameters given measurement sets. The 
next step in SARAS is development and implementation of such a code that 
operates on simulated data sets that include the non-ideal instrument behavior 
described in this work. 

We have designed SARAS with a 180° hybrid as a power splitter for the 
correlation receiver. The current implementation uses non-identical low-noise 
amplifiers. The analysis presented herein shows that using a matched pair of 
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amplifiers cancels certain contaminants. In the response to reference noise and 
receiver noise the imaginary components are canceled and, therefore, the net 
response of the system is now purely real. The real component is also altered 
in that this component of the response to reference noise is now ideal. The 
advantage of using matched amplifiers in the SARAS configuration is that the 
contamination in the real component of the calibrated spectrum, which con- 
tains the antenna temperature, is marginally reduced. However, the downside 
is that the calibration product in the form of the imaginary component of 
the response is lost. On balance, we have preferred to implement the system 
without using matched amplifiers because the additional contamination that 
this results in is expected to be only of a few mK amplitude. 

The system configuration could have used a quadrature hybrid that is a 90° 
power splitter in place of the 180° hybrid. Differencing the measurements in 
the two states of the cross-over switch would continue to cancel the correlator- 
based additive noise. In the case of a quadrature hybrid, the response of the 
system would manifest the antenna and reference temperatures in the imag- 
inary component instead of the real component. However, after calibration 
using the response to the calibration noise, the difference between antenna 
and reference temperatures would, as in the case of the 180° hybrid, appear 
in the real component. 

The primary advantage of using a quadrature hybrid is in the substantial 
reduction of reflection back to the input port when matched amplifiers arc 
used. This reduces the ripple in the response to antenna temperature because 
of the reduced reflection of this power back from the amplifiers to the antenna. 
The reference noise power, which dominates the ripple in the calibrated spec- 
trum in the case of a 180° hybrid, will continue to appear as the dominant 
ripple in the calibrated measurement set. Additionally, ripples owing to multi- 
path propagation of receiver noise will also continue to appear. The advantage 
of using a quadrature hybrid is that higher order ripples arising from multi- 
ple reflections between the antenna and amplifiers are reduced. The analysis 
presented above will be simplified, and higher order terms derived here may 
be neglected, if a quadrature hybrid along with a matched pair of low noise 
amplifiers replaces the power splitter in future. The alternate method of sup- 
pressing reflections from the low noise amplifier inputs, while retaining a 180° 
hybrid, is to implement each of the amplifier modules as a combination of a 
a quadrature hybrid followed by a pair of matched amplifiers followed by a 
second quadrature hybrid. 

The SARAS design presented herein has adopted a 180° hybrid followed 
by wide-band amplifiers without any quadrature hybrid. This choice was mo- 
tivated by practicalities in that the performance of available wideband 180° 
hybrids is superior to that of quadrature hybrids. In so far as we know, 180° 
hybrids can be made to be wideband and with smooth characteristics because 
they contain simple transformers; however, quadrature hybrids require a high 
order transfer function that will have correspondingly high order character- 
istics. If this is placed in the top of the receiver chain the modeling of the 
internal system will have more parameters to solve for. The real answer to 
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whether a quadrature hybrid or 180° hybrid ought to be used, and whether 
matched amplifiers are preferred in the SARAS configuration, will come from 
modeling using an analysis pipeline that derives likelihood estimates for sky 
model parameters given measurement sets. This may separately analyze the 
alternate configurations considering real-world models for the components to 
ascertain the structure for which the EoR parameters are better constrained 
in their likelihood estimates; as mentioned above this is part of future work in 
the SARAS development. 



7 Summary 

SARAS is a ground based experiment for precision measurements of the cos- 
mic radio background at low frequencies over an octave band from 87.5 to 
175 MHz. Measurement of a wideband sky spectrum at these long wavelengths 
with the exacting requirements demanded for detecting any faint signatures 
corresponding to rcdshiftcd 21-cm from the epoch of reionization motivated 
the development of a new design concept including a calibration strategy, and 
understanding the non-ideal behavior of the system to the extent that the 
science goals demand. 

The antenna in the SARAS configuration is a frequency-independent Fat- 
dipole over an absorber ground plane. We argue for a correlation spectrometer 
with the power splitter preceded by a cross-over switch, which acts as a differ- 
ential radiometer providing measurements of the difference between antenna 
temperature and noise temperature of an internal reference termination. The 
calibration noise is injected at the antenna via a directional coupler so that 
the interconnects and system behavior are unaltered during calibration. Ana- 
log signal transmission from the compact antenna electronics to a base station 
at a remote location is via an optic fiber link; this was another essential design 
feature that provided exceptional reverse isolation in the signal flow path and 
hence avoids unwanted coupling between the signals in the two arms of the 
correlation receiver. 

We provide a detailed analysis of the signal processing including inevitable 
effects related to multi-path propagation within the receiver electronics. We 
discuss coupling paths for the signal flow and the propagation of noise from the 
antenna, from the reference termination and from the low-noise amplifiers. We 
provide a formalism that describes the effects of multiple internal reflections 
between the impedance mismatches encountered by the signals, expanding the 
reflections as infinite series. We use measurements of reflection coefficients of 
sample devices, as well as of a frequency-independent Fat-dipole antenna, to 
derive example responses to antenna noise, reference termination and amplifier 
noise. These serve as a guide to the complexity necessitated of the modeling 
that derives constraints on parameters that describe sky spectral features. 

A key advantage of the SARAS configuration is that the measurement set 
is in the form of complex spectra with sky power appearing wholly in one 
component, the real component. Spectral structure arising from multi-path 
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propagation due to internal reflections appears in both components and, most 
importantly, shares parameters between the real and imaginary components. 
The system design also provides derived calibration products that are once 
again devoid of sky temperature and provide data that constrain the param- 
eters describing the spectral contaminants. 

This paper has described the system design and all of the rationale in the 
system design. The analysis of the signal flow and processing lays the founda- 
tion for the use of the system for measurements of the cosmic radio background 
spectrum. Varying the length of the interconnecting cable between the antenna 
and receiver varies the delay r in the above analysis and moves the responses 
arising from multi-path propagation into different Fourier components of the 
measured spectrum; thus SARAS measurements of the sky radio spectrum 
with differing lengths for I would be preferred for different science. SARAS is 
being deployed for absolute and precision measurements of the spectrum of the 
cosmic radio background in the 87.5 to 175 MHz band and for spectral features 
that may arise from hydrogen reionization in the redshift window 7.1-15.2. 
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